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ISSUE:  Dispersion  of  sediment  material 
from  a  disposal  site  must  often  be  assessed  by 
Corps  dredging  operation  elements.  Predic¬ 
tions  are  needed  in  cases  where  the  disposal 
site  capacity  or  retention  is  to  be  assessed  for 
some  purpose.  In  certain  special  cases, 
dredged  material  is  placed  in  an  underwater 
trench.  Whether  the  dredged  material  is  con¬ 
figured  as  a  mound  or  placed  in  a  trench,  spe¬ 
cial  procedures  are  required  to  predict  shear 
stress  on  the  material  and,  ultimately,  sedi¬ 
ment  transport. 

RESEARCH:  The  Dredging  Research  Pro¬ 
gram  (DRP)  includes  developing  guidance 
and  procedures  for  predicting  movements  of 
dredged  material  placed  in  open  water.  An 
important  part  of  DRP  research  is  to  check 
predictive  techniques  against  data,  as  was 
done  here. 


SUMMARY:  A  simple  numerical  model  is 
described  to  predict  velocity  profiles  in  the 
vertical  direction  and  bed  shear  stresses  at  var¬ 
ious  locations  along  a  mound  or  trench.  The 
computer  program  HPROFILE  operates  on  a 
personal  computer  and  is  easy  to  use.  The 
model  is  intended  to  allow  engineers  to  make 
calculations  as  part  of  an  assessment  as  to 
whether  dredged  material  will  remain  or  move 
from  its  placement  site. 

AVAILABILITY  OF  REPORT:  The  report 
is  available  through  the  Interlibrary  Loan  Ser¬ 
vice  from  the  U.S.  Army  Engineer  Waterways 
Experiment  Station  (WES)  Library,  telephone 
number  (601)  634-2355.  National  Technical 
Information  Service  (NTIS)  report  numbers 
may  be  requested  from  WES  Librarians. 

To  purchase  a  copy  of  the  report,  call  NTIS  at 
(703)  487-4780. 
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Preface 


This  study  was  conducted  at  the  Hydraulics  Laboratory  (HL)  for  the 
Coastal  Engineering  Research  Center  (CERC),  U.S.  Army  Engineer  Waterways 
Experiment  Station  (WES).  The  work  described  herein  was  authorized  as  part 
of  the  Dredging  Research  Program  (DRP)  by  the  Headquarters,  U.S.  Army 
Corps  of  Engineers  (HQUSACE),  and  was  performed  under  Work  Unit  32590, 
"Calculation  of  Boundary  Layer  Properties  (Cohesive  Sediments)." 

Messrs.  Glenn  R.  Drummond,  Rixie  J.  Hardy,  Vincent  Montante,  and  John  J. 
Parez  are  DRP  Technical  Monitors  from  HQUSACE.  Mr.  E.  Clark 
McNair,  Jr.,  CERC,  is  Program  Manager  (PM)  for  the  DRP,  and 
Dr.  Lyndell  Z.  Hales,  CERC,  is  Assistant  PM.  Dr.  Nicholas  C.  Kraus,  Senior 
Scientist,  CERC,  is  the  Manager  for  DRP  Technical  Area  1,  "Analysis  of 
Dredged  Material  Placed  in  Open  Waters." 

The  study  was  conducted  by  Mr.  R-  C.  Berger,  Estuarine  Simulation 
Branch,  Estuaries  Division  (ED),  HL,  in  conjunction  with  the  Work  Unit 
Principal  Investigator,  Mr.  A.  M.  Teeter,  Estuarine  Processes  Branch  (EPB), 
ED.  This  report  was  prepared  by  Messrs.  Berger,  Teeter,  and  W.  Pankow, 
EPB,  under  the  direct  supervision  of  Mr.  W.  H.  McAnally,  Jr.,  Chief,  ED,  and 
the  general  supervision  of  Messrs.  F.  A.  Herrmann,  Jr.,  Director,  HL,  and 
R.  A.  Sager,  Assistant  Director,  HL.  Ms.  L.  Pace,  under  a  contract  agreement 
with  Mississippi  College,  Qinton,  MS,  prepared  Appendix  B.  The  graphics 
software  used  in  the  personal  computer  program  (HGRAPH)  was  developed  by 
Mr.  D.  W.  Hyde,  Structural  Mechanics  Division,  Structures  Laboratory,  WES. 
Mr.  Steve  Knowles,  EPB,  collected  the  flume  data. 

At  the  time  of  publication  of  this  report.  Director  of  WES  was 
Dr.  Robert  W.  Whalin.  Commander  was  COL  Bruce  K.  Howard,  EN. 


For  further  information  on  this  report  or  on  the  Dredging  Research 
Program,  please  contact  Mr.  E.  Clark  McNair,  Program  Manager,  at 
(601)  634-2070,  or  Mr.  Allen  M.  Teeter,  Principal  Investigator,  at 
(601)  634-2820. 


Summary 


This  report  reviews  a  simplified  numerical  model  to  calculate  the  velocity 
profile  in  the  vertical  direction  for  changing  bed  slopes  subjected  to  currents. 
Model  predictions  were  compared  to  flume  data  for  flow  over  a  mound  and  a 
trench  to  determine  if  the  model  would  be  a  reasonable  planning  tool  for  cal¬ 
culations  involving  dredged  material  disposal  sites. 

The  flume  data  used  were  a  difficult  test  for  this  model,  but  the  overall 
profiles  compare  reasonably  well.  Therefore,  the  model  is  judged  reliable  for 
planning  purposes  involving  dredged  material  di^osal.  Although  the  model 
will  predict  the  vertical  velocity  profiles  in  low  riverine-  and  tide-dominant 
flows,  wave  action  is  not  included. 


1  Introduction 


Background 

The  cost  effectiveness  of  material  disposal  from  maintenance  dredging 
operations  depends  upon  the  distance  from  the  dredging  site  and  also  upon  the 
hydraulic  suitability  of  the  site,  the  type  of  material  being  deposited,  and 
available  equipment.  The  deposited  material  may  be  subjected  to  currents  that 
can  cause  it  to  be  eroded  and  transported  away  from  the  site.  The  severity  of 
this  condition  as  well  as  the  fate  of  the  material  need  to  be  addressed  in  the 
early  planning  phase.  The  precision  of  the  planning  estimates  depends  on  desk 
studies  typically  performed  on  a  personal  computer.  In  order  to  estimate  the 
amount  of  scour  and  transport  at  the  disposal  site,  the  bottom  shear  stresses 
must  be  developed;  therefore,  the  vertical  velocity  profile  within  and  adjacent 
to  the  site  is  of  interest.  For  uniform  bathymetry  conditions,  a  typical  logarith¬ 
mic  velocity  profile  suffices;  however,  a  more  general  method  is  needed  for 
variable  depths.  This  report  outlines  a  method  to  calculate  vertical  velocity 
profiles  over  variable  depths,  and  presents  comparisons  with  flume  data  to 
indicate  the  limitations  of  the  computer  code.  The  model,  HPROFILE,  pre¬ 
dicts  the  vertical  velocity  profile  at  various  depths  in  the  low  riverine-  and 
tide-dominant  flows,  but  does  not  include  wave  actions. 


Discussion 


In  a  uniform  flow,  it  is  customary  to  assume  the  vertical  velocity  profile  to 
be  logarithmic  as  shown  in  Equation  1: 


(1) 


where 


U  =  the  longitudinal  (horizontal)  velocity 


Chapter  1  Introduction 


1 


z  =  the  distance  above  the  bed 

u«  =  friction  velocity  (x/p)^^ 

T  =  shear  stress  on  the  bed 

p  =  fluid  density 

K  =  von  Karman’s  constant  (approximately  0.40) 

Zq  =  the  distance  above  the  bed  at  which  the  velocity  is  assumed  zero; 
the  roughness  length 

However,  when  the  bathymetry  changes,  the  flow  field  is  impacted  by  the 
modification  of  the  flow  boundary  layer.  Disposal  sites  may  indeed  have 
severe  depth  variations  such  as  trenches  or  formed  mounds,  and  it  is  important 
that  a  reasonable  velocity  profile  be  calculated  under  these  conditions.  The 
simplest  method  that  is  appropriate  is  van  Rijn‘s  method.^  The  basic  assump¬ 
tions  of  van  Rijn  are  that  the  vertical  velocity  profile  can  be  described  as  loga¬ 
rithmic  plus  a  perturbation  profile  due  to  the  bathymetric  effects;  and  that  the 
actual  profile  (with  bathymetric  effects)  and  the  equilibrium  profile 
(logarithmic)  are  identical  at  middepth.  The  manner  in  which  the  profile  is 
perturbed  to  account  for  the  changes  in  depth  is  empirical.  Details  of 
van  Rijn’s  method  are  given  in  Chapter  2. 

There  are  some  obvious  limitations  to  this  method.  The  first  is  that  it  is 
based  upon  empirical  data  (measurements)  and  may  compare  poorly  to  mea¬ 
surements  beyond  the  range  used  for  its  development.  Secondly,  it  is  stead)  - 
state.  This,  however,  is  probably  not  limiting  since  one  can  step  through  a 
tidal  cycle  with  fair  precision  as,  in  general,  the  temporal  terms  are  small. 

Van  Rijn,  on  the  other  hand,  was  able  to  get  reasonably  good  results  with  this 
simple  scheme.  In  fact,  it  compares  quite  well  with  a  much  more  sophisticated 
hydrodynamic  model  using  fc-t  turbulence  closure.  Some  of  van  Rijn’s  results 
are  shown  in  Figure  1  (where  this  simple  model  is  refened  to  as  the  profile 
model). 

Figure  1  indicates  that  the  method  is  capable  of  representing  flow  reversal 
due  to  an  expansion  of  a  mound.  However,  the  eddy  is  less  pronounced  and 
spread  over  a  larger  depth.  This  is  probably  due  to  the  limited  degree  of  the 
function  included  as  the  perturbation  profile.  Although  one  could  see  how 
assuming  that  the  middepth  velocity  matches  the  equilibrium  profile  could 
cause  this  effect,  it  is  a  remarkable  comparison  for  such  a  simple  model.  It 
also  demonstrates  the  relative  increase  in  lower  level  velocity  as  the  current 
flow  moves  up  a  slope  and  conversely  the  relative  reduction  of  the  deeper 
velocity  on  the  downward  slope. 


^  L.  C.  van  Rijn.  (1987).  "Mathematical  modelling  of  morphological  processes  in  the  case  of 
suspended  sediment  transport,"  Delft  Hydraulics  Communication  No.  382,  Delft  Hydraulics 
Laboratory,  Delft,  The  Netherlands. 


Chapter  1  Introduction 


Flume  Data 


Data  sets  with  bed  roughness  outside  the  range  originally  used  by  van  Rijn 
were  collected  in  a  Hydraulics  Laboratory  test  flume  at  the  U.S.  Army  Engi¬ 
neer  Waterways  Experiment  Station.  There  are  two  data  sets,  one  representing 
a  mound  and  one  a  trench.  The  flume  was  34  m  long  by  approximately  2  m 
wide,  and  the  testing  took  place  in  the  approximate  center  of  the  flume.  The 
mound  was  2.15  m  long  with  a  rise  of  0.051  m  and  1V:3.5H  side  slopes,  and 
the  trench  was  1.56  m  long  and  0.048  m  deep  with  1V;6.5H  side  slopes.  The 
mound  was  constructed  with  wood  and  the  trench  was  troweled  concrete  with 
surfaces  being  rather  smooth.  The  remainder  of  the  flume  surface  was  brushed 
concrete  texture. 

The  flume  data  collected  consisted  of  bottom  shear  stresses  and  cunent 
velocity  measured  over  the  vertical  as  well  as  horizontal  components.  Fig¬ 
ures  2  and  3  indicate  the  measured  bottom  shear  stress  in  the  flume  as  well  as 
the  bathymetry  of  each  configuration.  Figure  2,  the  trench  test,  indicates  that 
the  stresses  are  slightly  higher  downstream,  probably  due  to  the  increased 
bottom  velocity  as  flow  moved  up  the  slope.  Figure  3,  the  mound  test,  is 
more  difficult  to  interpret:  on  the  upstream  slope  (,ne  dip  in  stress)  prior  to 
the  mound  crest  (location  2.41  m),  the  flow  appears  to  be  at  a  stagnation  point; 
however,  this  is  not  the  case.  The  slope  is  1V:3.5H;  however,  the  plot  has  a 
distorted  scale  that  makes  the  slope  appear  much  steeper  than  it  really  is. 
Downstream  of  the  mound,  the  bottom  velocity  is  actually  lower  as  is  shown 
by  the  shear  stress.  The  spike  (location  0.098  m)  is  within  the  middle  of  a 
flow  reversal  close  to  the  bottom.  The  velocity  profile  data  collected  in  the 
flume  will  be  presented  with  a  comparison  to  numerical  model  data. 


Flume  and  Model  Comparisons 

The  numerical  model  (HPROFILE)  was  run  for  comparison  to  the  flume 
data  with  the  following  parameters: 

a.  Discharge: 

Q  =  0.104  m^/sec/m 

h.  Distance  above  bed  for  zero  velocity: 

2q  =  1  X  10'^  m 

Calculations  were  made  at  21  locations  for  the  mound  tests  and  22  locations 
for  the  trench  tests.  The  average  distance  between  stations  was  about  half  the 
depth.  The  model  requires  numerical  integration  from  station  to  station  in 
order  to  calculate  the  surface  velocity.  The  technique  used  is  a  Runge-Kutta 
fourth-order  scheme  so  that  there  is  some  stability  limit  on  how  large  the 
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distance  between  stations  can  be.  Details  of  the  model  and  its  application  are 
presented  in  Chapter  3.  An  analysis  of  the  stability  limit,  presented  in  Appen¬ 
dix  A,  demonstrates  that,  for  practical  purposes,  the  limit  is  sufficiently  ';rge 
that  it  should  not  interfere  with  normal  operations.  In  general,  the  distance 
between  stations  should  be  less  than  twice  the  depth,  and  the  size  of  adjoining 
cells  should  not  vary  too  quickly.  The  larger  of  the  two  adjoining  cells  should 
be  less  than  twice  the  size  of  the  smaller  since  gradual  transitions  are  more 
precise. 

Velocity  profiles  at  several  flume  locations  are  shown  in  Figures  4-19.  Fig¬ 
ures  4-10  are  the  mound  test  comparisons  and  Figures  11-19  are  the  trench 
tests.  The  figures  progress  in  the  downstream  direction.  It  should  be  noted 
that  the  model  is  presently  programmed  to  handle  only  positive  velocities  (i.e., 
velocity  in  the  direction  of  positive  Jt).  In  order  to  compare  the  flume  data,  the 
model  output  direction  was  reversed.  The  upstream  profile  is  assumed  to  be 
an  equilibrium  profile,  and  Zq  was  adjusted  until  the  upstream  reproduction 
was  good.  Figures  4-6  show  the  change  in  profile  as  the  flow  moved  upslope. 
The  model  seemed  to  react  more  quickly  than  the  flume  data  suggest,  becom¬ 
ing  uniform  sooner  than  the  flume  results.  Station  1.98  m  (Figure  7) 
improved.  However,  as  the  flow  moved  down  the  back  slope  (Figures  8-10), 
the  model  was  unable  to  pick  up  the  near-bottom  flow  reversal  (Figure  9).  It 
is  not  surprising  that  this  is  the  case  since  the  separation  was  so  near  the  bed 
and  since  the  perturbation  profile  includes  only  two  additional  functions.  At 
the  downstream  end  (station  -0.30,  Figure  10)  the  profile  compares  more 
closely,  but  it  is  evident  that  the  model  was  reacting  more  quickly  to  bed 
changes  than  was  the  flume.  The  flume  profile  is  more  uniform  in  the  surface 
layers  than  the  model  is. 

The  tests  conducted  provided  severe  conditions  for  the  model  to  reproduce. 
The  velocity  gradients  were  steep  and  confined  to  a  small  region,  and  the  per¬ 
turbation  profile  could  not  reproduce  these  conditions  very  well.  When  the 
profile  is  more  gradual,  the  model  does  relatively  well.  The  transitions  with 
the  bathymetric  changes  appear  somewhat  fast  in  the  model.  Any  correction 
for  this  would  need  to  be  made  via  the  empirical  constants  used  to  predict  the 
change  in  surface  velocity.  However,  not  enough  data  are  available  at  this 
time  to  be  sure  of  this  conclusion  nor  to  make  the  appropriate  modification. 
Overall  the  comparisons  are  sufficiently  precise  for  the  model  to  function  as  a 
planning  tool. 
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2  Details  of  Van  Rijn’s 
Formulation 


The  basic  premise  of  the  method  is  that  a  nonuniform  flow  profile  can  be 
described  as  a  linear  combination  of  a  logarithmic  profile  representing  the  law 
of  the  wall  and  a  perturbation  profile  representing  the  influence  of  pressure 
gradients: 


Uiz)  =  ^f/;,  In 


/  \ 
2 


+  Uv 


1  -  A  In 


[2n 


(2) 


where 

U(2)  =  velocity  at  distance  z  above  the  bed 
A,  t  =  constants  for  a  particular  longitudinal  location 
Uf^  =  velocity  at  the  water  surface 
h  =  depth 

n  =  (2  -  2o)/(/l  -  2o) 

This  yields  three  unknown  quantities  (A,  t,  requiring  three  equations  for 
solution. 

The  calculation  of  A  uses  conservation  of  mass: 

h 

Q  =  b  j  U(z)  dz  (3) 

^0 

where  b  is  the  width. 
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Substituting  Equation  2  into  3  and  integrating  provides  the  relationship  to 
calculate  A,  given  Q: 


The  calculation  for  parameter  t  includes  an  assumption  about  the  profile 
shape.  According  to  the  velocity  proGles  measured  in  a  trench  by  van  Rijn, 
the  middepth  velocity  is  approximately  equal  to  the  equilibrium  middepth 
velocity  given  by  the  logarithmic  profile: 


Finally,  the  quantity  U|^  is  calculated  by  the  relationship: 

^  -a  -a 

dx  h  ^  h 


where  jc  is  the  longitudinal  coordinate  and  is  the  surface  velocity  for 
equilibrium  flow  given  by 


The  calculations  of  parameters  and  a2  are  given  by  empirical  relation¬ 
ships  derived  from  flume  experiments: 


6 
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a,  =  0.28  +  0.11  tanh  6  (£L  -  O.is) 

'  [  iJ 


a,  =  0.235  0.065  tanh  17  _  -  0.035 

^  dx 
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3  The  Model  Code 


The  calculation  of  A  and  t  is  straightforward,  but  the  integration  of  the 
differential  equation  requires  some  additional  effort: 


=  Oi 


h,e 


“2  — 


(10) 


This  numerical  integration  was  canied  out  using  the  Runge-Kutta  fourth-order 
scheme.  In  this  method,  multiple  steps  of  the  relationship  were  performed: 


=  KVh) 


(11) 


where  /(f//,)  is  defined  as  the  right-hand  side  of  Equation  10.  The  steps  are  as 
follows.  Let  Vq  =  3t  the  initial  position: 


a.  Step  1:  Vj  =  Vq  +  1  Ax/(Vo) 

4 

b.  Step  2:  V2-V^*l  Axf{V^) 


(12) 


(13) 


c.  Step  3:  V3  =  V2  *  1  Ax/(V2) 


(14) 
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d.  Step  4:  V4  =  V3  + 


(15) 


UtSx„  -  Ax)  =  V4 


(16) 


This  method  is  very  precise  and  stable  as  long  as  the  length  step  Ax  is  not 
large. 

In  this  model  the  upstream  node  is  assumed  to  have  a  uniform  flow 
velocity;  therefore,  the  model  should  be  started  from  a  location  where  this  is 
true. 
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4  Summary 


In  this  report  a  simplifled  numerical  model  to  calculate  tlie  vertical  velocity 
proHle  for  changing  bed  slopes  was  reviewed.  Model  predictions  were  com¬ 
pared  to  flume  data  for  flow  over  a  mound  and  a  trench.  The  model  is  a  rea¬ 
sonable  planning  tool  for  calculations  involving  dredging  or  dredged  material 
disposal  sites. 

The  flume  data  used  were  a  difficult  test  for  this  model.  However,  the 
overall  profiles  compare  reasonably  well;  therefore,  the  model  is  judged  reli¬ 
able  for  planning  purposes.  Although  the  model  will  predict  the  vertical  veloc¬ 
ity  profiles  in  riverine  and  tidal  flows,  wave  action  is  not  included. 

In  order  to  effectively  use  this  program,  the  following  information  is 
required  by  the  equations  in  Chapters  2  and  3  and  should  be  collected  within 
and  adjacent  to  the  proposed  dredging  or  dredged  material  disposal  site(s): 

a.  Depth  and  tide  range. 

b.  Velocity:  surface,  middepth,  and  bottom. 

c.  Bathymetric  conditions. 
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HEIGHT,  Z,  m  HEIGHT. 


0.)  mis 

t— I  VERTICAL  VELOCITY 
0  2  nvs 

I - 1  longitudinal  velocity 


COMPUTED  longitudinal  FLOW  VELOCITY  OF  K-E  MODEL 

COMPUTED  LONGITUDINAL  FLOW  VELOCITY  OF  PROFILE  MODEL 

COMPUTED  VERTICAL  FLOW  VELOCITY  OF  PROFILE  MODEL 

EQUIVALENT  ROUGHNESS  .  0.02  m 

VON  KARMAN  CONSTANT  .  0  35 

DISCHARGE  .  0  08  m‘/s 

MEASURED 


Figure  1 .  Comparison  of  the  present  model,  a  k-t  model,  and  measured  data  sets  according  to 
van  Rijn 
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-0.1 


DISTANCE  ABOVE  BED.  M  «§  DISTANCE  ABOVE  BED,  M 


VELOCIPi'  M/S 


re  4.  Vertical  velocity  profile  for  mound  test  at  3.35  m 


Figure  5.  Vertical  velocity  profile  for  mound  test  at  2.68  m 
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e  6.  Vertical  velocity  profile  for  mound  test  at  2.41  m 


Figure  7.  Vertical  velocity  profile  for  mound  test  at  1.98  m 
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e  10.  Vertical  velocity  profile  for  mound  test  at  -0.30  m 


Figure  1 1 .  Vertical  velocity  profile  for  trench  test  at  3.35  m 
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re  12.  Vertical  velocity  profile  for  trench  test  at  2.59  m 


Figure  13.  Vertical  velocity  profile  for  trench  test  at  2.43  m 
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re  14.  Vertical  velocity  profile  for  trench  test  at  2.26  m 


Figure  15.  Vertical  velocity  profile  for  trench  test  at  1.74  m 
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16.  Vertical  velocity  profile  for  trench  test  at  0.61  m 


Figure  17.  Vertical  velocity  profile  for  trench  test  at  0.52  m 
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Figure  19.  Vertical  velocity  profile  for  trench  test  at  0.0  m 


Appendix  A 
Model  Stability 


Since  a  portion  of  the  model  involves  spatial  integration  to  calculate  the 
surface  velocity  f/^,  the  model  must  include  a  stability  limit  on  the  size  of  a 
spatial  step.  This  appendix  investigates  this  area  and  concludes  that  the  limit 
is  great  enough  that  reasonable  model  application  should  not  be  affected.  For 
simplicity,  consider  the  problem: 


(Al) 


where 


a2  =  parameter 
h  =  depth 

Solve  for 


(-^2 


U^(x)  =  U^io)  e 


(A2) 


The  numerical  representation  of  this  calculation  is  made  by  the  fourth-order 
Runge-Kutta  method.  This  multistep  calculation  is  made  as  follows  (where  for 
notational  simplicity  U  will  be  used  to  represent  and  V  represents 
intermediate  results): 


1/  rr  1  AJC  jj 


(A3) 
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The  subscript  n  stands  for  the  value  at  the  location  x  =  nl^  where  Ax  is  the 
step  length.  These  steps  may  be  combined  into  a  single  step  as; 


U, 


n+1 


(02  ^  -  4)(a2  ^  -  3)(a2  -  2)(a2  ^  - 1) 

24 


U, 


(A7) 


For  convenience  in  discussing  the  growth  of  round-off  errors,  assume  that 
Ug-  1.  It  is  assumed  that  the  round-off  error  enters  the  system  initially  with 
all  other  subsequent  calculations  being  exact.  The  "perturbed"  problem  will 
then  be  written  as: 


(02  —  -  4X02  —  -  3)(02  —  -  2)(o2  -  1) 

U*  ^  ^  ^  ^  ^  ru  (A8) 

n+1  24 


U* 

n 


and 


U*  =  1  -  e 
o 


(A9) 


where  e  is  the  initial  round-off  error. 

Subtracting  Equation  A8  from  A7,  defining  =  U^-  U*,  and  replacing 
the  multiplier  with  the  symbol  K  gives 

-  KE„  (AlO) 
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and 

£o  =  e 

(All) 

The  growth  of  the  error  can  be  traced  back  to  the  initial  state  since 
about  through  the  repeated  multiplication  by  K: 

it  comes 

E„=K\ 

(A12) 

The  initial  enor  should  decay  as  movement  occurs  away  from  the  origin; 
therefore: 

\K\  <  1 

(A13) 

which  in  turn  requires  that 

ct2—  <5 
^  h 

(A14) 

and  a2  has  maximum  and  minimum  values  of  0.30  and  0.17,  respectively. 

Selecting  the  most  restrictive  value  of  02  =  0.30  gives 

<  16.16 
h 

(A15) 

It  appears  that  relatively  large  spatial  steps  can  be  taken.  However,  the 
smaller  steps  would  be  more  accurate;  and  sirce  this  is  a  simple  one¬ 
dimensional  model,  high  resolution  is  possitie.  The  spatial  step  should  not  be 
a  severe  restriction. 

Appendix  A  Model  Stability 

A3 

Appendix  B 

HPROFILE  User’s  Guide 


Software  Installation 

Step  1.  Requirements  for  the  installation  of  HPROFILE; 

a.  IBM-PC  (AT  or  XT)  compatible  personal  computer  (PC),  with  either 
80286  or  80386  processor,  640K  of  memory,  MS/DOS  version  3.0  or 
later  operating  system,  VGA  or  EGA  color  monitor.  A  math  copro¬ 
cessor  is  recommended  for  faster  operation,  but  is  not  necessary. 

b.  Hard  disk. 

c.  Diskette  (floppy  disk),  5.25-  or  3.5-in.  with  the  following  files: 

HPROnLE.EXE 

DEMO.DAT 

Step  2.  Turn  on  the  computer  or  use  the  MS/DOS  command  to  change 
directory  (CD\)  to  the  system  prompt  as  shown; 

C:\>CD\  <RETURN> 

Step  3.  Using  the  DOS  command  MD  (Make  Directory),  create  a  new 
subdirectory  in  the  root  directory  for  the  HPROFILE  program  as  shown: 

C:\>MD\HPROFILE  <RETURN> 

Step  4.  Next,  type  CD\HPROFILE  <RETURN>  to  move  into  the  newly 
created  subdirectory.  The  prompt  will  now  be  C:\HPROFILE. 

Step  5.  Place  the  diskette  into  Drive  A  or  B.  Use  the  DOS  copy  com¬ 
mand  to  copy  the  required  files  onto  the  subdirectory  as  shown: 

C:\>HPROnLE  COPY  A:*.*  or  C:\>HPROFILE  COPY  B:*.* 
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Step  6.  The  message  "2  files  have  been  copied"  will  be  displayed  after  the 
files  have  been  copied.  To  verify,  use  the  DOS  command  to  look  in  the 
directory  by  typing  DIR  and  <RETURN>. 

Step  7.  HPROFILE  is  now  installed  in  the  computer.  Remove  the  diskette 
from  the  drive,  and  store  in  a  safe  location  away  from  the  computer.  The  rest 
of  this  guide  provides  additional  guidance  on  the  model  operation. 

Model  Operation 

Step  1.  To  execute  HPROFILE,  move  to  the  HPROFILE  subdirectory. 

Step  2.  HPROFILE  is  a  completely  menu-driven  system.  Type  HPROF, 
and  the  main  menu  will  appear  on  your  monitor  with  the  following  three 
options: 


Step  3.  To  make  a  selection,  press  the  corresponding  key,  which  is  in 
brackets.  Note:  Either  uppercase  or  lowercase  letters  can  be  used. 

Main  Menu  Option  Descriptions 

<R>un  model 

Before  the  model  can  be  executed,  the  data  file  must  be  prepared.  The  data 
file  should  be  an  ASCII  file  with  the  following  four  columns  of  data:  location 
of  nodes  (x),  distance  above  bed,  bed  slope  dhidx,  and  zero  velocity  bed 


B2 
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elevation.  The  data  file  should  be  on  the  hard  disk  of  your  computer  or  loaded 
on  a  floppy  disk  drive.  An  example  input  file  is  given  in  Table  Bl. 

Immediately  after  the  run  option  is  selected  in  the  main  menu,  the  follow¬ 
ing  model  parameter  menu  will  appear  with  the  default  values  for  the  data  set 
DEMO.DAT  entered: 


To  change  a  parameter  value,  enter  the  number  of  the  option;  the  cursor 
will  erase  the  old  value.  Enter  the  new  value  followed  by  <RETURN>. 
While  entering  parameters  please  make  note  of  the  following: 

a.  The  project  name  will  be  used  as  the  name  of  the  result  file  generated 
by  HPROFILE.  Make  sure  no  existing  file  has  this  name  because 
HPROFILE  will  destroy  the  file’s  contents  when  generating  the  result 
file. 

b.  The  project  name  and  input  file  name  must  be  less  than  12  characters 
long  including  .EXT 
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c.  The  units  listed  by  the  values  aie  the  only  units  HPROFILE  will  use. 

d.  There  may  be  no  more  than  25  nodes  in  one  execution. 

e.  Any  wrong  entry  can  be  changed  by  choosing  that  entry  a  second  time 
and  entering  the  correct  value. 

When  you  are  satisfied  with  the  model  parameters  entered,  choose  selection 
6  to  run  the  model.  The  message  "Executing  model"  will  flash  on  a  blue 
screen. 

Next  you  will  be  asked  if  you  wish  to  calculate  vertical  velocity  profiles. 
Enter  Y  for  yes,  and  the  prompt  "Node  =  "  will  appear  on  another  screen. 
Enter  the  number  of  the  node  the  vertical  velocity  profile  will  be  calculated  for 
followed  by  <RETURN>.  You  will  now  be  asked  for  the  number  of  points 
along  the  vertical  for  which  velocities  should  be  calculated.  Enter  the  number 
and  press  <RETURN>.  This  number  should  be  less  than  25.  The  profile  will 
be  echoed  on  the  screen  along  with  the  shear  stress  at  the  node  selected. 

When  you  are  ready  to  continue,  press  C  <RETURN>. 

You  will  be  asked  again  if  you  wish  to  calculate  a  vertical  velocity  profile. 
Enter  N  for  no  when  you  have  calculated  all  the  profiles  you  wish  to  for  the 
current  parameter  set.  The  main  menu  will  reappear  on  the  screen  and  the 
result  file  will  be  complete  and  will  be  saved  under  the  project  name.  An 
example  results  file  is  given  in  Table  B2. 

<G>enerate  graphics 

Two  types  of  graphs  are  available  with  HPROFILE;  x-y  plots  of  vertical 
velocity  profiles  and  x-y  plots  of  the  shear  stress  at  each  node  along  the  range 
of  the  bed.  When  the  graphics  selection  is  made,  the  following  menu  will  be 
displayed  on  the  screen: 


When  a  graph  type  has  been  selected,  the  current  project  name  will  appear 
on  the  screen  and  you  can  choose  to  accept  the  project  or  change  to  a  new 
project.  To  change  projects,  enter  C  and  the  cursor  will  erase  the  old  project 
name.  Enter  the  new  project  name  and  press  <RETURN>.  The  project  name 
must  be  the  name  of  a  result  file  generated  by  HPROFILE.  When  the  project 
name  is  acceptable,  press  A  to  accept  the  project  and  continue. 

To  generate  a  vertical  velocity  profile  plot,  press  V  <RETURN>.  The 
project  screen  will  appear.  Enter  the  correct  project  name  or  accept  the  current 
project.  The  next  screen  will  allow  you  to  enter  parameters  for  the  plot  and  is 
called  the  "Vertical  Velocity  Profile  Plot  Menu." 


VERTICAL  VELOCITY  PROFILE  PLOT  MENU 


(1)  No.  O'f  pro-files  on  plot  O 


<2)  Plot  title  Vertical  Velocity  Plot 


(3)  X  axis  ladel  Velocity  m/sec 


(4)  Y  axis  label  Height  m 


(5)  All  values  correct...  Draw  plot 


Enter  choice; 


Choice  1  allows  you  to  specify  how  many  profiles  will  appear  on  the  plot. 
You  may  have  up  to  five  profiles  on  one  plot.  The  data  for  each  profile  are 
generated  during  execution.  Only  profiles  for  nodes  specified  individually  at 
the  "Calculate  vertical  velocity  profile"  prompt  are  generated;  therefore,  only 
profiles  for  those  nodes  may  be  plotted.  After  you  have  entered  the  number  of 
profiles  to  appear  on  the  plot,  another  menu  will  appear  that  allows  you  to 
enter  the  node  numbers  for  each  profile  on  the  plot  and  the  description  that 
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will  appear  in  the  legend  (must  be  25  characters  or  less).  When  you  have 
entered  the  values  for  each  profile,  you  will  be  asked  if  these  values  are  cor¬ 
rect.  If  they  are  not  correct,  choose  the  appropriate  response,  and  you  will  be 
asked  if  you  want  to  change  a  node  or  legend  value.  Type  N  <RETURN>  to 
change  a  node  value  or  L  <RETURN>  to  change  a  legend  value.  Next  you 
will  be  asked  which  node  or  legend  value  you  wish  to  change.  Enter  the  cor¬ 
rect  number  and  you  will  be  allowed  to  make  the  change.  Be  careful  not  to 
accidently  press  the  <ESC>  key  or  any  function  keys  as  they  may  distort  the 
appearance  of  the  screen.  When  all  values  are  correct,  choose  the  appropriate 
selection,  and  the  "Vertical  Velocity  Profile  Plot  Menu"  will  return  to  the 
screen. 

Choice  2  allows  you  to  enter  the  plot  title.  It  may  be  up  to  40  characters 
long.  Choices  3  and  4  allow  you  to  label  the  x-axis  and  the  y-axis.  Regard¬ 
less  of  the  labels  entered  in  this  menu,  velocity  is  plotted  in  units  of  meters  per 
second  along  the  x-axis,  and  the  vertical  level  of  the  velocity  is  plotted  in 
meters  along  the  y-axis.  Any  wrong  entries  can  be  corrected  by  selecting  the 
entry  a  second  time  and  entering  the  correct  value.  Choose  5  to  draw  the  plot. 
A  message  will  appear  on  the  screen  asking  where  the  plot  should  be  drawn; 
choose  D  and  the  plot  will  appear  on  your  monitor.  If  for  some  reason 
HPROFILE  cannot  find  the  vertical  velocity  profile  section  of  the  result  file  of 
the  specified  project  for  a  particular  node  requested,  the  message  "Some  data 
not  available"  will  appear  on  the  monitor  and  HPROFILE  will  return  to  the 
graphics  menu.  Check  the  result  file  to  verify  which  nodes  profiles  were 
generated  for.  These  are  the  only  nodes  that  can  appear  on  the  vertical 
velocity  profile  plots.  Figure  B1  shows  an  example  vertical  velocity  profile 
plot. 

To  choose  a  shear  stress  plot,  press  S  <RETURN>.  The  project  selection 
will  appear,  and  after  you  have  selected  the  correct  project,  the  "Shear  Stress 
Plot  Menu"  will  appear. 


SHEAR  STRESS  PLOT  MENU 


(1)  plot  title  Shear  Strees  Plot 

<21  X  axis  label  Node  number 

<31  Y  axis  label  Shear  Stress  N/sg  m 

(4)  All  values  correct. . .Draw  plot 


Enter  choices 
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The  shear  stress  plot  di^lays  the  shear  stress  at  each  node  along  the  bed  on 
the  y-axis  against  the  node  number  along  the  x-axis.  To  change  the  plot  title 
or  axes  labels,  choose  the  appropriate  selection.  To  correct  an  entry,  select  the 
wrong  entry  and  enter  the  correct  value.  When  all  values  are  satisfactory, 
enter  4  to  view  the  plot.  Enter  D  in  response  to  the  prompt  for  the  plotting 
surface  and  the  plot  will  appear  on  your  monitor.  If  the  message  "Some  data 
not  available"  appears  on  the  screen,  check  the  project  result  file  to  see  if  it 
contains  the  proper  amount  of  data  for  each  node  in  the  range.  Figure  B2 
shows  a  sample  result  file. 

When  the  plots  (velocity  and  shear  stress)  appear  on  the  screen,  the  options 
"quit,"  "read,"  "zoom,"  and  "print"  also  appear.  To  find  out  the  y-value  for  a 
given  x-value,  press  R  or  move  to  the  "read"  option  using  the  arrow  keys  and 
press  <RETURN>.  Then  enter  the  x-value  of  the  noint  you  want  and  press 
<RETURN>.  The  corresponding  y-values  for  each  curve  will  be  printed  at 
the  bottom  of  the  screen.  Press  <RETURN>  and  the  original  four  options  will 
be  displayed. 

To  enlarge  a  specific  area  of  the  plot,  press  Z  or  move  to  the  "zoom" 
option  using  the  arrow  keys.  Press  <RETURN>  and  choose  "window"  and 
press  <R£TURN>  again.  Imagine  a  rectangle  encompassing  the  area  to 
enlarge  and  place  the  cursor  on  one  comer  of  the  rectangle;  press 
<RETURN>.  Move  the  cursor  to  the  opposite  comer  of  the  rectangle  and 
press  <RETURN>  again  (this  process  is  called  "rubber  banding").  The  screen 
will  be  redrawn  with  the  portion  of  the  plot  enclosed  in  the  rectangle  or 
"window"  enlarged  in  the  plotting  area.  You  may  use  zoom  again  to  enlarge  a 
portion  of  this  window  or  you  may  use  the  "zoom,  all"  option  to  return  to  the 
original  plot. 

The  "print"  option  executes  a  screen  dump  to  a  graphics  printer.  To  exit 
the  plot,  press  Q  and  the  program  will  return  to  the  HPROFILE  graphics 
menu. 

The  "Quit"  selection  in  the  graphics  menu  will  return  HPROFILE  to  the 
main  menu. 
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Table  B1 

EXAMPLE 

INPUT  DATA 

FOR  HPROFILE 

DISTANCE 

ZERO  VELOCITY 

X  LOCATION 

ABOVE  BED 

DH/DX 

BED  ELEVATION 

-3.3528 

0.317 

0.0 

.000001 

-2.98704 

0.317 

0.0 

.000001 

-2.68224 

0.317 

0.0 

.000001 

-2.52984 

0.317 

0.0 

.000001 

-2.48716 

0.3109 

-0.2857 

.000001 

-2.41402 

0.29001 

-0.2857 

.000001 

-2.3516 

0.27217 

-0.2857 

.000001 

-2.3165 

0.2662 

0.0 

.  00000 1 

-2.2000 

0.2662 

0.0 

.000001 

-1 .9812 

0.2662 

0.0 

.000001 

-1 .6764 

0.2662 

0.0 

. 00000 1 

- 1 . 3500 

0.2662 

0.0 

.000001 

-1 .0668 

0.2662 

0.0 

.000001 

-0.6096 

0.2662 

0.0 

.000001 

-0.3048 

0.2662 

0.0 

.000001 

-0. 1978 

0.2662 

0.0 

.OOOOOl 

-0. 1578 

0.27191 

0.2857 

.000001 

-0.09754 

0.28913 

0.2857 

.000001 

-0.02 

0.31129 

0.2857 

. OOOOO 1 

0.02 

0.317 

0.0 

.000001 

0.1524 

0.317 

0.0 

.000001 

0.3048 

0.317 

0.0 

.000001 

Table  B2 


SAMPLE  RESULT  FILE  HPROFILE 

Project  dew 

Input  ^ile  «  dew.dpt 

There  trt  22  1 

iKations 

Width  (!)  ‘  1.0 

Oischarqe  leu  i/wc)  *  .104 

NODE 

I  LOCATION 

DISTANCE 

SURFACE 

A 

T 

ZERO  BED 

SHEAR 

ABOVE  BED 

VELOCITY 

ELEVATION 

STRESS 

1 

-3.352800 

.317000 

.356196 

.078947 

.174050 

.000001 

.1147329 

2 

-2.987040 

.317000 

.356521 

.080573 

.174050 

.000001 

.1146161 

3 

-2.682240 

.317000 

.356739 

.081667 

.174050 

.000001 

.1145374 

4 

-2.529840 

.317000 

.356834 

.082140 

.174050 

.000001 

.1145034 

5 

-2.487160 

.310900 

.357039 

.048050 

.174228 

.000001 

.1194153 

6 

-2.414020 

.290010 

.358517 

-.075394 

.174870 

.000001 

.1376804 

7 

-2.351600 

.272170 

.360798 

-.190102 

.175462 

.000001 

.1553841 

8 

-2.316500 

.266200 

.362512 

-.226465 

.175671 

.000001 

.1615578 

9 

-2.200000 

.266200 

.367887 

-.195870 

.175671 

.000001 

.1595817 

10 

-1.981200 

.266200 

.376794 

-.147101 

.175671 

.000001 

.1563077 

11 

-1.676400 

.266200 

.386992 

-.094011 

.175671 

.000001 

.1525586 

12 

-1.350000 

.266200 

.395608 

-.051296 

.175671 

.000001 

.1493914 

13 

-1.066800 

.266200 

.401547 

-.022916 

.175671 

.000001 

.1472080 

14 

-.609600 

.266200 

.408826 

.010738 

.175671 

.000001 

.1445324 

15 

-.304800 

.266200 

.412459 

.027094 

.175671 

.000001 

.1431967 

16 

-.197800 

.266200 

.413550 

.031946 

.175671 

.000001 

.1427959 

17 

-.157800 

.271910 

.416869 

.084926 

.175471 

.000001 

.1351127 

18 

-.097540 

.289130 

.419546 

.202110 

.174898 

.000001 

.1163182 

19 

-.020000 

.311290 

.420155 

.323238 

.174217 

.000001 

.0963084 

20 

.020000 

.317000 

.418597 

.345491 

.174050 

.000001 

.0922566 

21 

.152400 

.317000 

.413715 

.327535 

.174050 

.000001 

.0940152 

22 

.304800 

.317000 

.408578 

.308181 

.174050 

.000001 

.0958655 

VERTICAL  VELOCITT  PROFILE 

NODE 

^  10 

NURBER  OF  POINTS  ALONG  THE  VERTICAL  ^  10 

ZERO  VELOCITT  BED  ELEVATION  < 

.000001 

1  LOCATION 

DISTANCE 

ETA 

VELOCITT 

ABOVE  BED 

1 

-1.98120 

.026621 

.10000 

.38608 

2 

-1.98120 

.053241 

.20000 

.40115 

3 

-1.98120 

.079861 

.30000 

.40467 

4 

-1.98120 

.106481 

.40000 

.40395 

5 

-1.98120 

.133100 

.50000 

.40116 

6 

-1.98120 

.159720 

.60000 

.39723 

7 

-1.98120 

.186340 

.70000 

.39262 

8 

-1.98120 

.212960 

.80000 

.38758 

9 

-1.98120 

.239580 

.90000 

.38227 

10 

-I. 98120 

.266200 

1.00000 

.37679 

SHEAR  STRESS  < 

.1563077  N/$q  ■ 
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